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We analyze both analytically and numerically the resonant four-wave mixing of two co-propagating 
single-photon wave packets. We present analytic expressions for the two-photon wave function and 
show that soliton-type quantum solutions exist which display a shape-preserving oscillatory exchange 
of excitations between the modes. Potential applications including quantum information processing 
are discussed. 
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I. INTRODUCTION 

The cancellation of resonant linear absorption and 
refraction via electromagnetically induced transparency 
(EIT) 0, has led to a range of new possibilities in 
non-linear optics. One important application is optical 
frequency mixing close to atomic resonances which allows 
the use of the enhanced nonlinear interaction without suf- 
fering from linear absorption and refraction. It has been 
predicted that EIT could lead to a new regime of nonlin- 
ear optics on the level of few light quanta 0, 0, IE IE ■ 
Several schemes for resonant nonlinear processes have 
been proposed and analyzed, both theoretically and ex- 
perimentally lEl- A particularly interesting system is 
resonant four-wave mixing using atoms with a double- 
lambda configuration |E m> 0| • Efficient frequency 
conversion, generation of squeezing 1121 as well as the 
possibility of mirrorless oscillations jl3l | with extremely 
low thresholds and narrow linewidth have been predicted 
|14| and in part experimentally observed [ISj]. 

Most theoretical and experimental studies of resonant 
nonlinear processes have been done for classical fields. 
For one-dimensional set-ups where common co-moving 
frames exist, full analytical solutions have been derived 
for the interaction of classical pulses in the adiabatic limit 
|16|. In view of the potential for an efficient nonlinear 
interaction on the level of few photons, however, a full 
quantum theoretical analysis of these systems is neces- 
sary. In addition, in order to take into account finite size 
effects which become increasingly important in the few- 
photon regime, and to analyze the potential for quantum 
information processing, such a quantum analysis has to 
go beyond the usual linearization approaches. Here very 
little work has been done, the few exceptions being the 
integrable models of resonantly enhanced Kerr interac- 
tion and photon blockade [E 0| ■ 

In a previous paper we have shown numerically that 
if single-mode fields are considered, it is possible to use 
an atomic vapor in a four-wave mixing double-lambda 
configuration to obtain full conversion from two input 
fields into two generated fields within a few centimeters 
of interaction length, even if the input fields only consist 
of single light quanta. The treatment was fully quantum 



but being restricted to a single- mode analysis, propa- 
gation effects of wave packets were not considered. In 
the present paper we extend this study by undertaking a 
multimode analysis in one spatial dimension and for co- 
propagating fields. In order to keep the problem tractable 
we restrict ourselves to the important special case of two 
single-photon wave packets as inputs. For this case we 
obtain simple analytic solutions and compare them to nu- 
merical simulations. We show that soliton-like quantum 
solutions exist that display an oscillatory exchange of ex- 
citation between the two input and the two generated 
fields. Finally we discuss briefly possible applications of 
single-photon four-wave mixing to quantum information 
processing and entanglement generation. 



II. SYSTEM AND EFFECTIVE FIELD 
EQUATIONS 

The situation we consider is resonantly enhanced 
four-wave mixing in the modified double-lambda system 
shown in Fig. ^ Note that a five-level atomic system 
is used instead of the original four-level system put for- 
ward in [E ^3 ■ This is due to the fact that in the four- 
level system the finite detuning A is associated with an 
AC Stark effect, which leads to intensity dependent dy- 
namical phase shifts of the fields. These phase shifts are 
of minor consequence in the case where the fields are 
counter-propagating ^E|> but for co-propagating fields 
have a detrimental influence leading to impaired phase 
matching and inefficient frequency conversion. As shown 
in 1^ these phase shifts are eliminated in the flve-level 
scheme when the relative sign of the dipole moments of 
the |4) |2),|1) transitions is opposite to that of the 
|3) — > |2),|1) transitions. 

The two fields with carrier frequencies vqi and h'n2 and 
slowly varying amplitudes fli and ^2 are initially excited 
and form the pump fields. The other fields of carrier 
frequencies i^ei and and slowly varying amplitudes 
El and E2 are generated during the interaction process 
and are assumed to be initially zero. Hi and Ei are taken 
to be exactly on resonance, i.e. vqi = UJ52, = W51, 
while the other two fields are detuned by an amount A, 
i.e. vn2 = ^ii - A = W31 -I- A and ue2 ^ ^i2 - A = 
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FIG. 1: Four- wave mixing in a modified double- A system with 
sgn[ci42/d4i] ~ — sgn[d32/ci3i], with dij being the dipole mo- 
ment of the jj) — |j) transition. 



ujr^2 + A. A finite detuning A, large compared to tlie 
Rabi frequencies, the Doppler broadening and the decay 
rates from the excited states, is necessary to maximize 
the ratio of nonhnear gain to hnear absorption. Decay 
from the two lower levels is considered to be negligible 
and all fields have the same propagation direction. 

Because of energy conservation there is overall four- 
photon resonance, i.e. i^qi + uq2 = ^ei + ve2- It can 
be shown that the contributions of the resonant transi- 
tions to the linear refractive index vanish if the fields are 
pairwise in two-photon resonance. Phase matching will 
thus favor two-photon resonance and we assume that this 
condition is fulfilled for the carrier frequencies of the four 
pulses, i.e. vei - mi = va2 - ve2 = '^21- 

Extending the analysis of .7j to a multimode descrip- 
tion, the interaction can be described by the effective 
adiabatic Hamiltonian [23 



^ hgc r n\alEiE2 + E\E\nih2 

n-mt - /dZ ^-j-^ — A+A ■ 



n\Vti + E\Ei 



Here CIi{z),Vl\{z) etc. denote dimensionless, slowly- 
varying (both in time and space) positive and negative- 
frequency components of the corresponding electric fields 
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L is the quantization length and k > —v/c is the wave- 
vector component in z-direction relative to p/c. The four 
fields are assumed to have either sufficiently different car- 
rier frequencies or different polarizations. Thus operators 
corresponding to different fields commute. The commu- 
tator between positive and negative frequency compo- 
nents can be approximated by a spatial delta function 



k{z 



Sij S{z - 



(4) 



We furthermore assume in |^ that all four transitions 
give rise to the same coupling g — Nd^uoi/ {2hceQ) — 
SNX^j/Stt, where N is the atomic number density, A 
the typical wavelength of the fields and 7 the typical ra- 
diative decay rate. It should be noted that numerator 
and denominator of Eq. ^ commute and thus there is 
no ambiguity with respect to the ordering of the terms. 

The structure of the denominator results from the sat- 
uration of the two-photon transition |1) — |2), whose co- 
herence lifetime is taken to be infinite. If a finite decay 
rate 70 of the |1) — |2) coherence is taken into account, a 
term proportional to 770 has to be added in the denom- 
inator. 

The non-polynomial character of the interaction 
Hamiltonian causes the nonlinear coupling to behave un- 
usually: As shown in (3], the interaction increases with 
decreasing pump field strength, making effective nonlin- 
ear frequency conversion possible even for single pho- 
tons. In the derivation of the effective Hamiltonian in 
adiabatic conditions were assumed. This limits the 
applicability of (^) in the multimode case to sufRciently 
long pulses. A discussion of non-adiabatic corrections 
and their effect on the propagation of the pulses is, how- 
ever, outside of the scope of the present paper and will 
be discussed elsewhere. 

As shown in the Appendix the slowly varying ampli- 
tudes of the electric field obey 

dt Ej {z, t) = -cd, Ej (z, t)+'- [Hint, Ej (z, t)] (5) 

and similarly for f2j. Thus from we arrive at the 
following equations of motion 

{dt + cd:}jEi = i^^ck(^\hlElElE2-n\^hlhlh2^k 
dt + cd^^(li = iKcA(^El^Clinin2- eI^IEiEiE2^A 
dt + cd^^E2 = -jKcA£;Jfii02 
(dt + cd;,'jQ2 = ~iKckn\EiE2 (6) 

where k = g/A and A = {n{ni + eIEi)-'^ . 

These equations admit four independent constants of 
motion: 

{dt + cd,)in{ni + ElEi) - 
idt + cd,){nlCi2 + ElE2) - 

{dt + cd,){Cl\Cli~CllQ2) ^ 



{dt + cd,){n\nlEiE2 + nin2EiEi) = o, (7) 

which represent the quantum analogs of the Manley- 
Rowe relations and additionally the conservation of the 
relative phase between the fields. 

If we assume that the input fields consist of two, single- 
photon wave packets in Qi and fl2, then it is clear, due to 
the constants of motion |71l , that the state of the system 
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can be represented at all times by 

(8) 

where \nk mk' Pk" Qk'") denotes n photons in the fcth 
mode of ili, m photons in the fc'th mode of CI2, and 
so on. 



III. FIELD INTENSITIES 

Since for the case of a single-photon input the expec- 
tation values of all fields vanish at all times, all relevant 
information about the state of the system is given by the 
mean intensities of the fields 

Mt)\n]iz)n,{z)\ip{t)), {ip{t)\El{z)E,{zMt)), (9) 

and the two-photon wave functions 

Mz,z',t) = {o\niiz)n2iz')\if{t)) 

oc 

^ ^2^^kz/L^2^^k'z'/L^^^^^^~^^^^Q^ 

k,k' — — 00 

^PEiz,z',t) = {0\Eiiz)E2{z')\m) 

00 

^2mkz/L ^2mk'z'/L^^^^^^y(^^^^ 



k,k' — — Qo 



ijjfi{z, z' ,t) and ipE{z,z',t) represent the amplitude of 
finding the Cli (Ei) photon at position z and simulta- 
neously the ft2 {E2) photon at position z' . The -(/j's are 
2D Fourier transforms from the fc-space representations 
^kk' and rjkk' into a real-space representation. 

We first discuss the dynamics of the mean intensities of 
the fields. Due to the constants of motion it is sufficient 
to calculate say (f2|f2i). 

(dt+ cd,yip\h\h^\ip) = {ip\'- [i?,f)Ir!i] \^) 

= iKc{{h\VLih\h\EiE2 - (lih2^ElVL\n[^k) 

= iKc {eIeICiiCi2 - n\nlEiE2) (12) 

where we have dropped the common spatial coordinate 
z and used -ff|0) = 0, A\(p) = \ip) as well as the fact 
that three annihilation operators acting on \(p) give zero. 
Similarly we find 



(dt + ca,) {ElElnin2 - n\n\EiE2 



= 2iKc{Vl{VllVLiVL2 - e\eIEiE2) (13) 
{dt + cd^^ {(i\VlIq.i(12 - E\eIEiE2) 

^ 2iKc{E\E\VLi(l2 - Cl{nlEiE2)- (14) 
Consequently the differential equation we must solve is 



We take the input to be two independent single-photon 
wave packets in fii and Q2 with the same spatial en- 
velope fo{z) and vacuum in the other two fields, which 
corresponds to a separable initial state of the form 



®|0)^;i®|0) 



E2' 



The are Fourier-transforms of fQ{z) 



z/L 



Thus 



{nl{z,t)ni{z,t))in = {nl{z,t)n2{z,t))i 

= ipo{z~ ct) 



(16) 



(17) 



(18) 



where i^o{z) — fo{z)- With these initial conditions one 



finds 



= {(11^2) = i^o{z - ct) cos^{kz), (19) 



i.e. a sinusoidal exchange of excitation between the two 
pump and the two generated fields. Complete conversion 
is achieved at z = tt/k. It is worthwhile noting that 
as shown in 0, and contrary to the classical dynamics, 
complete conversion can only be achieved in the quantum 
case for initial Fock states with one or two photons in the 
two pump modes. 



IV. DYNAMICS OF THE TWO-PHOTON WAVE 
FUNCTION; QUANTUM SOLITONS 

We proceed by calculating the two-photon wave func- 
tions. The calculation of 'ipa,{z, z', t) can be split into two 
distinct cases, depending on whether z = z' or z ^ z' . 
Let us first assume z = z' . 



dt+cdA{Q\nin2\^) - (o| 



-iKc{0 I hih2{(l\VLlEiE2 + VLiQ.2E{Et)^^\ 
-iKc{0\EiE2W (20) 



where we have again used the results that H\0) = 0, 
A\ip) — \ip) and that three annihilation operators acting 
on \ip) give zero. Operating on {0\EiE2\4') we find 

idt+cd,){0\EiE2\^) = -ik(0|M2|¥'>- (21) 

Thus, to determine the wave function, we need to solve 
the differential equation 



(dt + cd,'^ {n\ni) = -4K^c^{dt + cd,yCl\ni). (15) {dt+cd,'^ i^n{z,z,t) ^ -K^c^i^niz,z,t). (22) 
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For an input consisting of two independent single-photon 
wave packets with the same spatial shape, as considered 
above, the initial two-photon wave function reads 



Im (z,z,t) ] 



i>n{z, z', t)in = f{z - ct)f{z' - ct). 
We also have at the entrance of the medium 
[dt ^ cdAi^n{z,z,t) = -iKc{{)\EiE2\f) 

SO that the solution to (|22|) is given by 



z=0 



(23) 



(24) 



^^{z,z,t) = '4)o{z - ct)cos{nz), (25) 

with t^oiz) = f'^{z). Thus the two-photon wave function 
at equal spatial points propagates through the medium 
modulated by a factor of cos(kz). We see that after one 
full conversion cycle z — t:/k the phase of the wave func- 
tion has changed sign. This agrees with a numerical sim- 
ulation of the quantum problem, the results of which are 
shown in Fig. |21 




Z / K " 



FIG. 3: Generated fields in 45-mode simulation. Both initial 
wave packets taken to be Gaussian. 



Re[V|/^(z,z,t)] 




t/(cK) 



Z / K" 



FIG. 2; 45-mode simulation of wave function propagation in 
a co-moving frame. Both initial wave packets taken to be 
Gaussian. 



We can use a similar procedure to find ^E{z,z,t) — 
{Q\EiE2W), the wave function of the generated fields. 
We obtain 

ipsiz, z, t) = —i ijjQ{z — ct) sin(Kz). (26) 

The evolution of this field is shown in Fig. |3| Note that 
it is tt/2 out of phase with the drive field wave functions, 
as expected. 

For full knowledge of the state of the system, we also 
need to find ?/'(z,z',t) where z ^ z' . From we see 
that operators at different points in space commute, so 
that 



dt+ cd, + cdA{Q\vii{z)vi2{z')\>p) 



(27) 



Expanding the Hamiltonian, normal ordering the expres- 
sion and again noting the form of we eventually arrive 
at 



dt + cd,+cdA (0|f2i(z)f22(z')|(^) =0. 



(28) 



Thus the wave function of the system in the case where 
z ^ z' is simply given by the corresponding input expres- 
sion 



^P{z,z',t) = fiz-ct)fiz' ~ct), z^z'. 



(29) 



(0| i/,f2i(z) tl2{z')^^x{z) H,Vl2{z') \^) 



It is evident that the initial two-photon wave function at 
different coordinates propagates undisturbed throughout 
the medium, and that there is a discontinuous change 
in the behavior when moving away from the line z = 
z' . This is clearly seen in Fig. 01 which shows the two 
dimensional wave function ip^i^z, z' ,t) at time r = tt/k. 
This is the time required for one full conversion from 
the pump fields to the generated fields and back again. 
Essentially we have at this time 

■il:n{z,z',T) ^ -ipfi{z,z')in z = z' (30) 
■00(2, z',r) = '(/'o(2;,2;')in z^z' (31) 

This behavior can easily be understood. The two- 
photon wave function represents the joint probability 
density that two photons are located at positions z and 
z' . As both photons propagate at the same speed they 
will never meet at the same point in space unless they did 
initially. Since the nonlinear interaction in |^ is local, 
two photons at different spatial points do not interact. 
Thus only if z = z' is there a conversion from the pump 
fields into the generated fields due to the local nonlinear- 
ity. 

The simple sign change in the wave function for z — z' 
hints at the possibility of using the system as a phase 
gate for quantum computation, as was mentioned in the 
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Re[\|/^(z,z',T)] 



Z'/ 




'ipn{z,z') ^ = (l)o{z) S{z - z') 



(33) 



Z / K" 



FIG. 4: Numerical 45-mode calculation of ipsn{z, z' ,t), shown 
after one full cycle of energy transferal from pump to gener- 
ated fields and back again. The sign flip for 2; = 2' is evident. 



single- mode case considered in 7]. One chooses the 
length of the nonlinear medium such that exactly one 
full conversion cycle can occur. If the entire wave func- 
tion changes sign, we would have a system that behaves 
as a true phase gate: if only one of the two inputs is 
populated, the pulse exits the medium unchanged, but if 
both are present, a phase shift of tt occurs. 

However, due to the multimode nature of this problem, 
it is ip{z,z\t) for all z and z' that is a true reflection 
of the system rather than ^{z,z,t) and there is no sign 
change for z ^ z' . To see whether we can still use the 
system as an approximate phase gate, we consider the 
behavior of the state vector coefficients in ((H)). Returning 
briefly to a picture with 2N + 1 discrete modes, by using 
Eq. I|38() of the Appendix one can show that (|30|l and 
(|5T|l imply 



2N 

■mn 



emn(O). (32) 



This gives a good picture of how the different modes have 
mixed among themselves, and shows that in general a 
simple sign change of all the coefficients cannot exist. 
Only if we can enforce that at least temporarily only a 
single effective mode of the two pump fields is excited, 
e.g. by using a resonator set-up is it possible to use the 
system as a phase gate. 

We see however from (|32() that after one full conversion 
cycle the initially factorized state H16|l evolves into an 
entangled state between all modes of the fields Vli and 0,2- 
Thus the nonlinear interaction generates entanglement. 

On the other hand if the initial state is a two-photon 
wave packet in an entangled state, such that the initial 
two-photon wavefunction has only a contribution for z = 



only diagonal components of the two-photon wave func- 
tion will ever be nonzero. According to H25() and H26() 
they undergo sinusoidal oscillations 

ipn{z,t) = (l)o{z - ct) cos{kz), (34) 
i/)£;(z,t) = — i0o('Z ~ ct) sin(Kz). (35) 

The superposition of pump and generated fields 

$(z,t) = cos{Kz)'>pfi{z,t) + ism{Kz) ipE{z,t) 

= Mz-ct) (36) 

propagates in a form-stable manner and represents a 
quantum soliton solution. The two-photon wavefunction 
$(z,t) of the quantum soliton corresponds to a quasi- 
particle excitation 

*t(^,t)= (37) 
cos{Kz)il\{z , t){}l{z, t) — i sin{Kz)E\{z, t)E\{z, t). 



V. SUMMARY 

In the present paper we have presented a full quan- 
tum theoretical treatment of resonant forward four-wave 
mixing for pulses. For this we have used an effective 
Hamiltonian derived from the interaction of the fields 
with an ensemble of atoms in a double-lambda configura- 
tion using an adiabatic approximation 0. We were par- 
ticularly interested in the few-photon regime since here 
quantum effects dominate the dynamics and because of 
potential applications in quantum information process- 
ing. Thus we have restricted ourselves to the important 
special case of an input consisting of two single-photon 
wave packets. For this case we were able to analytically 
solve the propagation equations for the field intensities 
and two-photon wave functions which contain all rele- 
vant information about the quantum state. 

We found that there is an oscillatory energy exchange 
between the two pump and generated fields with 100% 
conversion at periodic intervals of interaction. This result 
is characteristic for a few-photon Fock-state input; for a 
coherent input complete conversion can only be achieved 
asymptotically for very large input power. 

We have also shown that after even multiples of 
the conversion length the two-photon wave function 
■0(2, z' , t) regains its initial form, while after odd mul- 
tiples there is a sign flip for z = z' . 

If the two input wave packets are not independent but 
in a highly entangled state, the two-photon wave function 
can be made zero outside of the diagonal. It was shown 
that such a pair of input wave packets form a formstable 
soliton-like quantum solution which is a superposition of 
pump and generated fields with oscillating coefficients. 
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The process of resonant four-wave mixing was shown to 
generate large entanglement between the modes forming 
the two, single-photon wave packets. Furthermore the 
nonlinear interaction strength is large enough to gener- 
ate a controlled phase shift of a single photon by the 
presence of another one. Thus, if the number of relevant 
modes is at least temporarily restricted by some external 
means like a resonator, the system could have interesting 
applications as a photonic phase gate. 
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Appendix 

To simplify the transition from a single-mode descrip- 
tion to a multimode one, we will first consider one single- 
mode quantum field a, and then generalize the result 
to our four-field system. Suppose that the single-mode 
Hamiltonian governing the evolution of a is given by 
H = fkoo o^a + Hintia, a^). To go over to the multimode 
description we consider an interaction region of length 
L, divided into 2N + 1 cells, and consider a discrete set 
of modes around the carrier frequency of the field, i.e. 
kn = ko + 2mT/L, —N < n < N. We now define local- 
ized field operators (denoted by a tilde) via 



The multimode Hamiltonian can now be written as 



ai 



ak 



N 

E 

k=-N 



ak exp 



N 



2iV-|- 1 



2N - 



^ ai exp 



1 ^ 



~2Trikl 
2N+1 



(38) 
(39) 



where the cik are annihilation operators for mode k. a/ is 
related to the field strength in cell I. 

These localized field operators have the commutation 
relations 



ai,dl 



{2N + l)5u 



(40) 



k I 

hujQ 



2N - 



(41) 



I II 
where u)q is the carrier frequency and 

N 



E 

k=-N 



2nkc 



{2N + iyL 



exp 



2Trik{l - I') 



2N+1 



(42) 



Commutation relations yield the following Heisenberg 
equation of motion 



d ^ 

—ai = -uJoai 
at 



i{2N +l)^uju,ai, 
I' 



^[Hint,a] (43) 



Now, as we let iV ^ oo we find 

IL/{2N+1) -> z (44) 
o; d(z,t) (45) 
~ii2N+l)Y,^wai' -> ~cd,aiz,t) (46) 

[a(z,t),a^z',t)] -> LS{z-z'). (47) 
Introducing slowly time- varying amplitudes we obtain 



{dt + cd,)a{z,t) = -[Hin,{t),a{z,t)] 



(48) 



Thus the multimode equations of motion look exactly like 
the single-mode equations of motion with the exception 
that a cdz term has been added and the fields now have 
a spatial dependence. 

Returning to the four-wave mixing situation described 
in this paper, we see that the interaction Hamiltonian 
given by |^ is of the form shown above, with the sum- 
mation replaced by an integral in the continuum limit. 
The multimode annihilation and creation operators E(z) 
and rt{z) defined in and are analogous to the lo- 
calized field operators a{z,t) defined above, except for 
the factor of y/L inserted to ensure correct commutation 
relations regardless of the quantization length. Thus the 
equation of motion (jS)) follows from H48|l above. 
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